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["t I . The mixed states are important in quantum optics since they frequently appear in 



the decoherence problems. When one of the components of the system is prepared in 
the mixed state and the evolution operator of this system is not available, one cannot 
Qh' deduce the density matrix. We present analytical approach to accurately solve this 

S3 , problem. The approach can be applied on the condition that the Schrodinger's 

^ . ' equation of the system is solvable with any arbitrary initial state. In deriving the 

, ^! 

solution we exploit the fact that any mixed state can be expressed in terms of a 

^ , phase state. The approach is illustrated by deriving the density matrix of a single- 
in 

OO ' mode heat environment interacting asymmetrically with two qubits. Our results are 
^ ! 

■ in good agreement with the available results in the literature. This approach opens 

. new perspectives for treating complicated systems and may impact other applications 
in the quantum theory. 

^ ! PACS numbers: 03.65.Ud, 03.67.-a, 42.50.Dv 
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I. INTRODUCTION 

In quantum theory, the state of the system has two categories; namely pure state and 
mixed state. A pure state implies perfect knowledge of the system. For the mixed case. 



Electronic address: 



erorany@hotmail.com, Tel:006-166206010, Fax: 0060386579404 



2 



we do not have enough information to specify the state of the system completely and hence 
cannot form its wavefunction. In this situation, we can only describe the system via the 
density matrix pf. Mixed state frequently appears in the decoherence problems, e.g., 
We can obtain the mixed state, which is associated with any state, by considering its phase 
to be totally randomized. The best example of the mixed state is the thermal field, which 
represents an electromagnetic radiation emitted by a source at temperature T. The exami- 
nation of a quantum beam with a thermal noise is an important topic from both theoretical 
and practical points of view, e.g. {2]. 

In the quantum information theory, a considerable attention has been paid to the en- 
tanglement of the bipartite and the multipartite systems in which one of the subsystems 
exists initially in the thermal equilibrium 3|-|5|. The conclusion of all these studies is that 
it is possible for the thermal field, which is a highly chaotic field, to induce entanglement 
between qubits. The previous studies have been limited to the systems whose exact form of 
the evolution operators is obtainable js, 5|. For the other systems the solution is difficult or 
even impossible. For some few systems of the latter, one has to use the numerical methods 
to solve the master equation of the system 4[], but the results cannot be totally trusted. In 
this paper we develop, for the first time, a simple analytical approach solving this problem 
exactly. This approach can provide the dynamical density matrix of the unitary quantum 
system whose one of its components is initially in the mixed state such as the thermal field. 
It works only when the Schrodinger's equation of the system for any arbitrary initial state is 
solvable. Furthermore, the approach can be applied as an alternative technique for finding 
the solution of the systems whose evolution operators exist but they are very complicated. 
For instance, the evolution of the mixed field with the multi-level atoms, e.g. three-level 
atom and four-level atom, etc. 

In section II we describe the approach in details. In section III we give an application 
for the approach by solving specific problem. The example, which we have considered, is 
that the problem of the two atoms in the cavity interact asymmetrically with the thermal 
field. The reason behind this choice is that this system is a subject of the current research 



in relation to entanglement 



3|-|5|. Thus, we can easily validate the approach by comparing 



our results with the available results in the literature. 
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II. DESCRIPTION OF THE APPROACH 

In this section, we describe how one can deduce the density matrix of the unitary quantum 
system whose one of its components is initially prepared in the thermal state. Before going 



into details, let us briefly state some properties o 



a diagonal expansion in terms of Fock states as jl|: 

OO „ 

Pf = = - / \<^))W)\d4>. (1) 



the thermal field. The thermal field has 



n=0 







where p{n) is the photon-number distribution of the thermal field and |-2(0)) is the associated 
phase state having the form: 

OO OO 

\zm = VpRexp(m0)|n) = Y^Cnmn), (2) 

n=0 n=0 

where 

= \/pMexp(m0), p{n) = -^j-^-^ (3) 

and n is the mean-photon number of the thermal field having the form n = (e''"'/'^*'^— and 
kh is the Boltzmann's constant. It is obvious that n increases by increasing the temperature 
T. 

Now we are in a position to explain the approach. Assume that we have a bipartite system 
whose Hamiltonian is Hdb, where d and b stand for the field and the other party, respec- 
tively. These two components are initially prepared in the states \ipd(yO)) and lipbi^))- The 
Schrodinger's equation of this system is solvable and hence we can obtain the wavefunction 
as iV'(t)). From the basic concepts of the quantum mechanics we have: 

|^(t,0)) = ?7(t, 0)1^,(0)) (g) 1^,(0)), (4) 

where U(t, 0) is the evolution operator of the system regardless if we can deduce its explicit 
form or not. Under these assumptions we can solve this system accurately if the party d is 
initially prepared in the thermal field pf instead of {ipdi^))- In this case the initial state of 
the system reads: 

TT 

m = Pf^\Mo)){Mo)\ = \j \z{4>)){zm^m^)){umd(t>- (5) 
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Under the evolution operator U(t,0), the system in the initial state ([5]) evolves as: 

TT 

p{t) = f/(t,o)p(o)c/t(t,o) = (6) 



where [^/'^(t, 0)) = U{t,0)\z{(j))) (^{ipbiO))- Our target is to deduce the wavefunction 
\'^z{t,4'))- If the explicit form of U{t,0) is available, the solution is straightforward. If 
it is not, we can solve the Schrodinger's equation for the initial condition 1^(0)) and \tpb{0)). 
We have to emphasize if the system is solvable for the arbitrary state |V'rf(0)), it will be 
automatically solvable for the state \z{(f))). This is based on the fact that any state of the 
field is just a linear combination of the Fock states weighted by a specific distribution. As 
soon as we derive \ipz{t,4')) we substitute it into ([6]) and carry out the integration over the 
phase (f) to obtain the requested density matrix. From the above description, the theme of 
the approach we transform the problem of obtaining the dynamical density matrix to that 
of finding the wavefunction of the system when the field is initially prepared in the state 

\zm. 

In spite of the simplicity of the approach, it is efficient and able to provide the exact 
solutions for some non-analytically solvable problems so far. In the following section we 
show this fact by developing the analytical solution for one of those problems, which has 
been already numerically treated. 



III. EXAMPLE: TWO-QUBIT PROBLEM 

In this section, we apply the approach described in the preceding section to obtain the 
density matrix of the two-qubit problem. Precisely, we deduce the dynamical density matrix 
of the system of the single-mode thermal field interacting simultaneously with the two two- 
level atoms in the cavity in the non-identical fashion. The numerical technique has been 



exploited for solving this system {4], however, here we present the analytical solution. It is 
worthwhile mentioning that the two-identical-atom version of this system has been already 
de„.onst.ated, e.g., ffl. Fo. the latte. the dynamical den..ty tnattix of the system ha. 
been derived by the Kraus representation [7]. That cannot be applied to the non-identical- 
atom system. Also, we verify the validity of the approach by comparing some of our results 
to those in the literature in relation to entanglement. We proceed, under the rotating wave 
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approximation the Hamiltonian describing the system takes the form 4-|6|: 

^ = Hq + Hi, 

(7) 

2 

Ho = ua^a + + Hj = J2 Aj(aa+ + a^a'), 

where and Hj are the free and interaction parts of the Hamiltonian, and are the 
Pauh spin operators of the jth atom; a (a^) is the annihilation (creation) operator denoting 
the cavity mode, u and Ua are the frequencies of the cavity mode and the atomic systems 
(we consider that the two atoms have the same frequency) and Xj is the coupling constant 
between the jth atom and the field. Based on the relation between Ai and A2 we have 
two cases; namely asymmetric case (Ai 7^ A2) and symmetric case (Ai = A2). We should 
stress that the asymmetric case is closer to experiment than the symmetric one Finally, 
throughout the investigation we assume that coa = ou. 

Now, our aim is to deduce the explicit form of the density matrix p{t) of the Hamiltonian 
d?]) when the field is initially prepared in the thermal field ([T]) and the atoms are initially in 
the following mixed state: 

Pa = sin^6'cos^?9|ei,e2)(e2,ei| + sin^ 6'sin^ ^?|5fi, 5f2) (5-2, S-il + cos^ 6'|ei, 5f2)(5'2, ei|, (8) 

where \ej) and \gj) stand for the excited and the ground atomic states of the jth atom, 
respectively. The variables 6 and i) are phases, which can be specified to give different 
types of the initial atomic states. The subscript a denotes the atomic system. We should 
stress that the suggested approach is applicable for any initial atomic state not only (|8]). 
Noteworthily, for Ai 7^ A2 we cannot derive the explicit form of the evolution operator of the 
system U{t, 0), however, we can solve the Schrodinger's equation |[9|]. The latter is the main 
demand to apply the approach. The whole initial state of the system can be easily written 
in the form ([5]). Under the evolution operator U{t,0), this initial state evolves as: 

p(t) = f/(t,0)/)(0)f/t(t,0) 

^ ^ 
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where 

|vl/,(t,0)) = f/(t,O)|vI/,(O,0)), J = 1,2,3, 

|^i(0, 0)) = |ei, 62, ^(0)), 1^2(0, 0)) = \gu 92, ^(0)), 1*3(0, 0)) = |ei, ^2, 2:(0)). 



(10) 



Our main concern is to derive the wave functions |\l/j(t, 0)). These can be easily obtained 
by solving the Schrodinger's equation three times for the given initial conditions as 9|: 



*i(^,0))= ECn(0) Xi'\t,n) \ei,e2,n)+X^"{t,n) |ei,^2,n + l) 

n=0 



+ Xi'\t,n) \ g,,e2,n + l) +Xl'>{t,n) | ^1,^2^ + 2) 



0')/ 



J = 1,2,3 



where 



XI 



(1 
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X. 
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(1 
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(2 



X, 



(2 



X. 



(2 
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(2 



(3 



X. 
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X 



(3 



^(^> 


n) 
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n) 




^(i> 


n) 




\t, 


n) 




\t, 


n) 




\t, 


n) 




\t, 


n) 




\t, 


n) 






n) 




\t, 


n) 




\t, 


n) 






n) 





cos(tO+) _|_ + cos(fn„ ) 
^^■rl 2A„ r^n 2A„ ' 

:^%p{[4A?(n + 2) - /i-]^ - [4A?(n + 2) - 
^^^^{[4A2(n + 2) - - [4A2(n + 2) - 

2AiA2-v/(n+2)(n+l) r /j-O-M 
2AiA2A/n(n-l) 

yT— L'"^°V''"n-2; ~ >-^^V''"n-2;i 5 



in-2 

- jAi y/n 



[cOs(t^]+_2) - COs(til„_ 



[4A2(r2-l)+/i+_2; 



+ 1 sin(*"n-2) 



sin(tf2^_2) 

< I l—S- ^ _ I I # I I — I — II _ I I LS- /\ _ I II I I — I — II _ I 

2A„ 

-»A2v^ \\A\2(^ n I ,,+ 1 sin(m+_2) r/l\2,'^ i\ , ism(tfi;^_2) 



-[4Ai(n-l) + /i-_2] 



2A„-: 
M^-2 



[4A?(n-l)+/.+_2] 



-[4A2(n-l)+/i;_2] 



2A„_2 

iA2v^ 

2A. 



COs(tfi+_2 



/^n~2 
2A„_2 



COs(tl]„_2), 
sin(tn+_J 



r {K-i - 4A?(n + 1)]!^^^!^ + [4A?(n + 1) - /.ti 



I sin(ta„_i) 



i-^ {[A? - A2 + A„_i] cos(ta+_i) + {\l - A? + A„_i] cos(ta;;_i)} 
^|^[cos(tf]t,)-cos(tn;_i)], 
^ { + 4Ain] - + 4Ain] ' 



and 



(11) 



(12) 



An = V(A? + Ai)2(2n + 3)2 - 4(A? - \^^\n + 2)(n + 1), 



A? + A2 ± A, 



U± = ^V(A^ + A^)(2n + 3)±A„. 



(13) 



7 




(a) (b) 

FIG. 1: Atom-atom entanglement induced by the interaction with a thermal field having n = 1 
when the atoms are initially prepared in |ei,e2) (a) and |ei,g2) (b). 



The dynamical density matrix of the whole system can be obtained by substituting ffTTj) - 
(fT3|) into Q and then carrying out the integration with respect to 0. This treatment provides 
the exact solution of the problem. Moreover, it is easier and better than solving the Liouville 
equation numerically ^. To evaluate the density matrix of the two qubits, i.e. pa{t), one 
has to trace over the field variables the density matrix of the system as: 

(14) 

where 



J2 P{n) \x[^\t,n)\'^ sin' e cos^ ^ + \X['> {t,n)\'^ sin' e sin' ^ + \X['' {t,n)\'^ cos^ e 



(3), 



n=0 

oo 



B. 



^ P{n) \X^^\t,n)\'' sin' e cos'' ^ + \X^'> {1,71)1' sin' e sin' ^ + \Xl^> (t^n)]' cos' e 



9192 



(3), 



n=0 

oo 



Be,g2 = E Pin) Xi^^ (t, n)X^s''* (t, n) sin' 6 cos' 1^ 



(1)* 



J (15) 



n=0 



+ XP {t, n)Xf ^* (t, n) sin' 6 sin' ^ + Xf^ {t, n)Xf ^* {t, n) cos' 6 



(3)/ 



'(3)*/ 



B. 



B 



ei92' 



For the best of our knowledge this is the first time the explicit form of the density matrix 
of this system to be presented. This reflects the power of the approach. The density matrix 
f[T^ tends to that of the symmetric case [s^, which was obtained by the Kraus representation, 
by simply setting Ai = A2 = A. 

To verify the validity of the approach we comment on the entanglement of the two qubits 
controlled by ([1] 
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In this regard, we use the negative values of the partial transposition [10[, which is 



frequently used for such systems 



0,3 



, and is defined as: 

^ = -25^7-, (16) 
3 

where 7" are the negative eigenvalues of the partial transposition of Pait). The entanglement 
measure ranges between ^ = for separable qubits and ^ = 1 for maximally entangled qubits. 
For the density matrix (fT4l) there is only one eigenvalue of those of the partial transposition 
matrix, which can approach negative values, having the form: 



277 = i?e,e2 + B,,,, - y'(5e,e2 - B.^g,)^ + 4|5,,,,|2. (17) 

By means of (|T2l)-( IT71) we have obtained the results of for the same values of the inter- 
action parameters. Now it is therefore reasonable to make a comparison with the results of 
the general case, which has been numerically presented in jj]. In that article the analysis 
has been confined to such forms of the coupling constants Ai = 1 + 7 and A2 = 1 — 7, where 
7 (with < 7 < 1) is the relative difference between the two atomic couplings. This means 
that the strength of the interaction of one of the bipartites (atom-field) is increasing, while 
the other is decreasing simultaneously. As an example, in Figs. 1 we plot the quantity C, for 
the same values of the interaction parameters as those of Figs. 2 in [^J, which were given 



for the concurrence. It is obvious that our figures and those in {4] are identical even though 
they represent different measures. This clearly confirm the validity of the approach. Obvi- 
ously, the analytical treatments are in general better than the numerical ones as they can 
provide us some analytical facts about the system. For instance, the condition of involving 
the expression (fT7|l negative values is: 

= -Beiea-Bgiga - l-BgieaP < 0. (18) 

From ( fT2l) -( fT5|) . when n = and the two atoms are in \ei,g2), the inequality ( fTSi) can be 
simplified as: 

T = - ,,,^]%,, [cositnt,) - l][A?cos(^^]^l) + a^] < o. (19) 

It is evident that the expression f|T9|) gives T = sin^(tf2l^)/4 for the symmetric case. This 
means that the symmetric case (with n = 0) cannot generate any entanglement between the 
two atoms, which are initially in \ei,g2)- Nevertheless, the asymmetric case can generate 
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entanglement for certain values of the coupling constants, in particular, for those satisfying 
the inequality cos(tfil]^) < — Ag/ A^. Furthermore, for the atoms, which are initially in l^fi, (72) 
and n = 0, one can easily prove that Xf'{t, 0) = 0. Thus, the condition ([18]) is not fulfilled, 
i.e. entanglement cannot be established in this case. 

It is worth prompting that we have aimed by the above discussion to justify the validity of 
the approach not to repeat the study of the entanglement of the two qubits. So that we have 
selected few cases for the sake of comparison only. Nevertheless, from the treatment we have 
performed, which was not presented here for the sake of brevity, and the other treatments 
given in j3|-|5| we can conclude that a highly chaotic state in an infinite-dimensional Hilbert 
space can entangle two qubits depending on the type of their interaction with the field as 
well as the initial conditions of the system. The amounts of entanglement generated by the 
asymmetric case are much greater than those generated by the symmetric one. 

In conclusion, we have developed, for the first time, a simple analytical approach for 
solving unitary system when the initial field, as one of its components, is in the mixed state, 
e.g. thermal light. The approach is applicable only when the Schrodinger's equation of 
the system is solvable for any initial arbitrary state. We have verified the validity of the 
approach for some selective cases for the entanglement of two qubits. The approach enables 
us to check the results, which have been numerically obtained earlier 
we believe that our approach represents a powerful tool to solve such type of problems and 
may impact other applications in the quantum theory. 
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